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1. Introduction 

Accurate pressure broadened profiles of alkali resonance doublets perturbed by H 2 and 
He are of crucial importance for the modelling of atmospheres of late M, L and T type 
brown dwarfs and for generating their synthetic spectra in the region 0.6 - 1.1 /jm. 
The dominant lines are the Na I 589.0/589.6 nm and K I 766.5/769.9 nm doublets 
[H [2] which occur in the middle of a spectral region where background absorption is 
particularly small so that both the line centres and wings stand out. There can also 
be significant contributions from less abundant alkalis such as Li, Rb and Cs, and from 
subordinate doublets such as Na I 818.3/819.5 nm. 

The profiles of the strongly broadened alkali doublets have been the subject of 
several recent studies using models primarily designed to describe the non-Lorentzian 
far wings of the profiles. Burrows and Volobuyev [3] used a quasistatic unified Franck- 
Condon semiclassical model to study the Na I and K I doublets whereas a unified 
line shape semiclassical model has been used by Allard et al [I] for the Na I and K I 
doublets and, more recently, for the Rb I and Cs I doublets [5]. Alternatively, quantum 
mechanical single-channel models based upon bound-free and free-free transitions in the 
transient molecules formed during collision have been used to calculate the emission and 
absorption spectra for the wings of the Li I [HI [7] , Na I [SI IE] and K I [E] resonance lines 
in the presence of He perturbers. Despite these studies, highly accurate calculations of 
the central Lorentzian cores are still needed in order to estimate the effects of dust in 
brown dwarf atmospheres. 

We report results for the widths and shifts of the Lorentzian profiles over the 
temperature range 70 < T < 3000 K of the Li I, Na I and K I doublets and the 
Na I subordinate doublet broadened by helium perturbers. The calculations extend 
our previous study [9] of pressure broadening of the Na I doublets by He at laboratory 
temperatures up to 500 K and are based on a fully quantum-mechanical close-coupling 
description of the colliding atoms and the Baranger [10] theory of line shapes. As input 
for our calculations we have computed new ab initio potentials for the alkali-helium 
interaction. Preliminary results for the widths of two Na I lines for temperatures up to 
2000 K have been given in [TTJ and for the widths of the Na I and K I doublets in [12jjj. 

2. Interatomic potentials 

The adiabatic molecular potentials +1 Va(R) for the X*-He system (where X = Li, 
Na or K) have been obtained by using a three-body model in which the alkali X* is 
treated as a X + ion plus an active electron and the perturbing He atom is represented 
by a polarizable atomic core. Model potentials are used to represent the electron-atom 
and electron-atomic ion interactions and the basic methods adopted for obtaining these 
model potentials are discussed by [13j. These potentials reproduce not only the correct 
energies of the bound states of X* but also generate wavefunctions for both the bound 

§ Note that the results reported in 12J are for full half-widths and not half half-widths as stated 
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and scattering states of X+e~ and He+e - that contain the correct number of nodes. The 
model potentials also support unphysical bound states corresponding to the presence of 
closed shells in the X + ion and in He. This effect is taken into account in the calculation 
of the molecular potential by including the unphysical states in the atomic basis used 
in the diagonalization of the Hamiltonian for the three-body model. 

The most important recent change from the earlier approach is that the core-core 
interaction is itself calculated directly using a three-body model composed of X + and 
He + cores plus one electron. A more detailed discussion of this development is given for 
the Li- He case by [H]. Further numerical refinements have been incorporated into the 
computer program since then and new calculations carried out for the three alkali-helium 
systems considered in our study. 

The ns 2 Vs, np 2 Vs and np 2 Vn molecular potentials for Li-He (n = 2), Na-He 
(n = 3) and K-He (n = 4) are shown in figures [U [2] and [3] respectively. The Na-He 3d 
2 ^s,n,A potentials are shown in figure HI The minima in the potentials used and their 
positions are given in tables [T], [2] and [3] and compared with other data where they are 
available. New ab initio potentials for K-He have recently been obtained [27J but no 
details are given. The overall agreement with other work is very satisfactory and the 
three-body model can be expected to give good results at medium and large interatomic 
separations. However the positions of the minima in the np 2 Vn potentials occur at 
short range and are largely determined by the behaviour of the core-core interaction. 
For the Li-He and Na-He cases the agreement with theory and experiment for the np 
2 Vu states is generally good, but the result for the K-He 4p 2 Vn state requires some 
comment. Because the 3p virtual state supported by the K + + e~ model potential has 
a higher energy than that of He(ls 2 ), the core-core potential, K + -He, is unreliable at 
small separations. Methods for resolving this difficulty are still being investigated, but 
since the behaviour of line widths and shifts depends largely on the difference between 
upper and lower state potentials they are not particularly sensitive to this problem. 

3. Spectral line profile 

The width and shift of each spectral line have been calculated using the quantum- 
mechanical impact theory of Baranger [I10J for non-overlapping spectral lines in which 
the profile of each isolated line is a Lorentzian. The orbital, spin and total electronic 
angular momentum operators for the alkali atom are L, S and j = L + S respectively and 
radiation is emitted by the alkali atom as it undergoes a transition between initial and 
final states \riiLiSiji) and \rifLfSfjf), where n denotes the principal quantum number. 
The half half-width w and shift d at temperature T of the Lorentzian profile is given by 



[31j 




(1) 



where f{E) is the normalized Maxwellian perturber energy distribution 




(2) 
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iV is the perturber number density and 

x[5i,i> - (ji l ' Ji\S\ji I Ji)(jf I' Jf\S\j f I J f )*] (3) 

describes the effects of collisions on the two states forming the spectral line. Here M 
is the reduced mass of the emitter-perturber system, I and /' are quantum numbers 
corresponding to the relative emitter-perturber angular momentum L^j before and 
after the collision and J = L^j + j is the total angular momentum of the emitter- 
perturber system. The scattering matrix element in the coupled \j I J) representation 
is (f I' J\S\j I J) where we have suppressed the quantum numbers (n,L,S) for 
convenience and {j^/j} is the 6 — j symbol [32] . 

The scattering matrix elements in ([31) are determined from the asymptotic behaviour 
of the radial functions Gj t {R) for each scattering channel (j, I, J). These functions satisfy 
the coupled equations [3Tj 

[qjP - " ^2 ^ + k< j} G jlj"l"( R ) = V /lj'l>( R ) G j>l>,j"l"(R)> ( 4 ) 

where (j", I") labels the linearly independent solutions of (j3J) and the Born-Oppenheimer 
coupling terms have been neglected. The parameter 

fcj = 2M[E - (E^ + e,(oo))}/h 2 (5) 

is positive for open scattering channels and negative for closed channels. Here E is 
the total energy of the emitter-perturber system, E^ L ' S "> is the energy of the state of 
the separated atoms to which the molecular state dissociates adiabatically and the fine 
structure parameter Sj(R) has been assumed to have its asymptotic value e 3 -(oo). The 
interaction potential matrix elements Vjj dfl ,(R) are, for the ns 2 Sj level 

V jltjV {R) = 6 j)j >5 lil , 2 V E (R), (6) 
for the np 2 Pj levels, 

V jldV (R) = Sjj.Sut, 2 V U (R) + cV v fV s (R) - 2 V U (R)} (7) 
and, for the nd 2 Dj levelg[JJ] 

Vjtjn'iR) = hA,v 2 V±(R) + + (-l) l ' +l ]B jldn ^Vn(R) - 2 V A (R)] 

+C$U 2 Vx(R)- 2 V A (R)] (8) 

where the coefficients 

c fl]'v = E (-i) j '~ j c(j j i; -n n o)c(j f r ; -n n o) 

n 

xC(n \ j; O Q)C(n \ j'; O Q) (9) 

|| We correct here the phase factor given in [5] for the 2 T)j levels. 
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and 

Bjijui = J2 °( J J h -fi 0)C(J / -n n 0) 

xC(2 | j; 1 Q - 1 fi)(7(2 \ f; 1 - 1 fi) (10) 

are symmetric under (J, I) <->• (f,l'). Here Cijijuja, Qif2 2 fi 3 ) is the Clebsch-Gordan 
coefficient, f2 denotes the projection of an angular momentum onto the internuclear 
axis R and 2S+1 Va(R), where A = \Ql\, are the adiabatic molecular potentials. 

The equations (TjJ decouple into two sets of opposite parity (— 1) J±1 / 2 . Thus the 
ns 2 Sj, np 2 Pj and nd 2 Dj states give rise to one, three and five coupled differential 
equations respectively for each parity These equations were solved using a modified 
version [31 J of the R-Matrix method of Baluja et al [33] and the solutions fitted to 
free-field boundary conditions to extract the scattering matrix elements. 

4. Results and discussion 

Calculations have been completed for the Li I doublet 2p 2 P3/2,i/2 2s 2 Si/ 2 , the Na 
I doublet 3p 2 P 3 / 2 ,i/2 — > 3s 2 Si/ 2 , the K I doublet 4p 2 P 3 / 2 ,i/ 2 -> 4s 2 Si/ 2 and the 
Na I subordinate "doublet" 3d 2 D 3/2 -> 3p 2 P 1/2 and 3d 2 D 5/2i3/2 -»■ 3p 2 P 3/2 . The 
main computational issues associated with extending our previous calculations up to 
the higher temperatures (3000 K) of astrophysical interest are convergence of the sum 
over partial waves in ([3]) and of the integration (pQ) over perturber energies. The present 
calculations required about 500 partial waves for most of the energy nodes but up to 
1000 partial waves were needed at the highest energy nodes which corresponded to 
temperatures extending to 5000 K. 

The .R-matrix propagation was commenced using the arbitrary choice R = I 
at a distance -R m i n well within the classical turning point so that the solutions were 
independent of R m i n - Typically -R m i n ranged from 2ao for small I to 25a for large 
/. The distance R max at which the solutions were matched to the free-field solutions 
was typically 300a — 400a . This yielded S- matrix elements accurate to at least six 
significant figures and, after summing over partial waves, resulted in six figure accuracy 
for dtS(E) and four for ^sS(E). As the major part of the i?-matrix method is energy 
independent, equations (BJ are solved at each value of I for the entire set of energy nodes. 
The large number of partial waves needed might suggest that a semiclassical treatment 
would be more appropriate. However the present quantal calculations did not require 
an inordinate amount of computer time (typically a few hours) and are free of the 
uncertainties endemic in semiclassical calculations of determining the appropriate lower 
cutoff on impact parameters in order to exclude perturber paths entering the classically 
forbidden regions of the interatomic interaction [31]. 

The temperature dependence of the computed half half-widths and shifts are shown 
in figures [5] and [6] respectively. The widths and shifts of all lines were calculated at 10 
K intervals and in all cases were found to be smooth functions of temperature. The 
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widths are accurately represented to three significant figures by the power law form 

w(T) = aT b (11) 

where the fit parameters a and b are given in table HI Also shown are the values of b 
obtained by [35] using a semiclassical model which does not resolve the members of each 
doublet. Our results agree closely with this early calculation except for the anomalous 
4p 2 Pi/2 - 4s 2 Si/2 potassium transition where our value is significantly higher. For all 
lines the temperature dependence is stronger than the T 3 behaviour obtained from a 
classical treatment using a pure van der Waals V(R) = —Cq/R 6 interaction. 

Although the shifts have a more complex temperature dependence, they all have 
the general feature of increasing with temperature above T ~ 500 K apart from the 
4p 2 Pi/2 - 4s 2 Si/2 potassium transition which is again anomalous as its shift decreases 
with temperature. 

The effects of fine structure on the calculated widths and shifts arises from the 
6 — j symbols in ([3]) and the j-dependence of the S'-matrix due to the potential 
matrix elements Vjijqr(R) and the fine-structure splittings Ej(oo) appearing in (T4j). 
However, as the zero energy for the calculation of any level 25+1 Lj is set at the energy 
E = E( L ' S ^ + €j(oo) of that level and the splittings are negligible for all but the lowest 
energy nodes, the role of the actual splittings is relatively unimportant. This is evident 
in the results for the Na lines 3d 2 D 5 / 2 ,3/2 —* 3p 2 P3/2 where the widths and shifts are 
significantly different even though £3/2(00) = 65/2(00). 

The behaviour of the anomalous potassium transition warrants some comment. 
The calculation for the potassium 4p 2 Pi/2 level is different from that of all other levels 
considered in the present study in that there are closed channels at the lowest energy 
nodes and the wells of the molecular potentials are very shallow (depths E d < 10~ 4 au), 
smaller than the fine structure splitting e = 2.68 x 10~ 4 au of the level. Consequently 
the channels only become open at energies above e and therefore at energies above E d . 

We compare our calculated widths and shifts with measurement in table [5] and 
table [6] respectively. We have not included all the experimental width data for Na 
as the present calculations for the Na doublets closely reproduce our earlier results 
[9] for T < 500 K and a detailed comparison with the numerous existing theoretical 
and experimental studies was reported in that paper. In general the widths are in good 
agreement for all alkalis although only one measurement is for T > 700 K. The situation 
regarding the shifts is less clear. The predicted shift for the 3d 2 D 3 / 2 - 3p 2 Pi/2 Na I 
line agrees very closely with the measurement of [40] and the results for the K I doublet 
have the same sign and relative magnitude as the measurements. However the predicted 
shifts for the Li I and Na I doublets have the opposite sign to that of the measured shifts. 
The shifts are quite sensitive to the precise details of the potentials as they are produced 
by a balance between the effects of the long-range attractive potential and the short- 
range repulsive potential. In particular, for a given energy they are sensitive to where 
the repulsive wall is located. Consequently the disagreement with experiment for the Li 
I and Na I doublets suggests that the repulsive region of the Li potential may need to 



Alkali doublets perturbed by helium 



7 



be slightly suppressed and that of the Na potential enhanced. 
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Figure 2. Adiabatic np 2 V"e molecular potentials: Li-He ( ), Na-He ( — • — ) and 

K-He ( ). 
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Figure 3. Adiabatic np 2 Vn molecular potential: Li-He ( ), Na-He ( — • — ) and 

K-He (- -). Note that the Na-He and K-He potentials have been scaled up by factors 
of 5 and 20 respectively. 




Figure 4. Adiabatic potential curves for Na-He: 3d 2 Ve ( ), 3d 2 Vn ( — • — ) and 

3d 2 V A ( )■ 
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Table 1. Positions Rd (do) and minimum energies Vd (au) of the Li-He adiabatic 
molecular potentials. Also given are the Cq (au) coefficients for the asymptotic forms 
of the potentials. 



2s 


2 Vv 


2p 


2 Vx 


2p 


2 V n 


Rd 


v d 


Rd 


v d 


Rd 


v d 


11.73 


-7.005 x 10~ 6a 


16.48 


-2.596 x 10~ 6a 


3.39 


-4.470 x 10~ 3a 


11.553 


-6.787 x 10" 6b 


17.027 


-1.511 x 10~ 6c 


3.395 


-4.55 x 10~ 3d 


11.56 


-6.79 x 10- 6 ° 


18.0 


-1.4 x 10~ 6f 


3.44 


-5.65 x 10^3 g 


12.2 


-9.1 x 10~ 6k 


17.0 


-3.1 x 10~ 6h 


3.44 


-4.670 x 10" 3h 


12.0 


-4.6 x 10~ 6f 






3.44 


-4.07 x 10^3 1 


11.385 


-1.18 x HT 5h 






3.40 


-4.638 x 10~ 3j 


11.6 


-6.47 x 10~ 61 






3.60 


-3.14 x 10-3 k 










3.42 


-3.95 x 10" 3f 










3.37 


-4.647 x 10" 3 1 


C* 6 = 22.57 a 




C* 6 = 50.75 a 




C* 6 = 28.31 a 





a Theory, this work. 

b Theory QH OH [23] 

c Theory [J51[231[25] 

d Theory [HI EH [25] 

° Theory [TH] 

f Theory [29] 

s Theory [22] 

h Theory [TBI CIS] 

1 Theory [2BJ 

j Theory [2TJ 

k Theory [28] 

1 Experiment [2T] 
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Table 2. Positions Rd (ao) and minimum energies Vd (au) of the Na-He adiabatic 
molecular potentials. Also shown is the position R p (ao) and energy V p (au) of the 
peak in the 3d 2 Vs potential together with the Cq (au) coefficients for the asymptotic 
forms of the potentials. 



3s 2 ^ s 



3p 2 1/ E 



3p 2 14i 



R d V d 


Ra 


v d 


Rd 




v d 


12.14 -6.610 xl(T 6a 


18.77 
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-1.617 x 10~ 3a 


12.0 -1.3xl0~ 5b 
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43.72 a 




3d 2 y E 
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2 V U 




3d 
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Rd V d 
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3 3.08 -2.167 x 10~ 7a 

Rp Vp 

7.37 4.252 xl0" 3a 



-2.010 x 10" 
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-9.717 x 10" 



Ca = 218.7 £ 



C fi = 194.3' 



Ca = 121.2 : 



a Theory, this work. 

b Theory [IB] 

c Theory [22] 

d Theory [IB] 

e Theory [25] 

f Theory [2T] 

s Theory [T7j 

h Experiment [20] 
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Table 3. Positions Rd (do) and minimum energies Vd (au) of the K-He adiabatic 
molecular potentials. Also shown are the Cq (au) coefficients for the asymptotic forms 
of the potentials. 



4s 


2 Vv 


4p 2 l/ s 




4p 


2 V n 




V d 


R d V d 


Rd 




v d 


13.65 


-5.395 x 10" 6a 


20.52 -1.547 xl0~ 6a 


7.96 




-1.986 x 10~ 4a 


13.9 


-3.5 x l(T 6c 




5.29 




-2.19 x 10~ 3b 


11.5 


-2.85 x l(T 5c 




5.29 




-1.03 x 10~ 3c 








5.30 




-1.12 x 10" 3 d 








5.4 




-7.7 x 10~ 4e 


C* 6 = 42.24 a 




C 6 = 105.6 a 


0> = 


61.96 a 





a Theory, this work. 
b Theory [22] 
c Theory [21] 
d Theory [16] 
e Theory [30] 



Table 4. Parameters a and b for the fita) = aT b to the temperature dependence of 
the pressure broadened Lorentzian half half-widths (in units of 10 -21 MHz m 3 /atom) 
of various alkali doublets due to He perturbers. The results apply to the temperature 
range 70 - 3000 K except for the 4p 2 Pi/2 - 4s 2 Si/2 potassium transition where the 
fit applies to the range 500 - 3000 K. Numbers in parentheses denote the uncertainty 



in the last 


significant figure from the fitting 


I procedure. 






Element 


Transition 




A (nm) 


a 


b 


°LML 


Li 


2p 2 Pi/ 2 - 


2s 2 S 1/2 


670.78 


0.02533(5) 


0.3998(2) 


0.39 


Li 


2 P 2 P 3 /2 - 


2s 2 S 1/2 


670.79 


0.02461(5) 


0.4042(3) 




Na 


3 P 2 Pi/ 2 - 


3s 2 S 1/2 


589.59 


0.02011(6) 


0.4270(4) 


0.40 


Na 


3p 2 P 3 /2 - 


3s 2 S 1/2 


589.00 


0.02918(8) 


0.3866(4) 




Na 


3d 2 D 3/2 - 


3 P 2 Pi/ 2 


818.33 


0.06075(8) 


0.3799(2) 




Na 


3d 2 D 3/2 - 


3p 2 P 3 /2 


819.48 


0.06369(5) 


0.3737(1) 




Na 


3d 2 D 5/2 - 


3p 2 P 3 /2 


819.48 


0.05755(9) 


0.3820(2) 




K 


4p 2 P 1/2 - 


4s 2 S 1/2 


769.90 


0.01385(8) 


0.4886(8) 


0.39 


K 


4 P 2 P 3/2 " 


4s 2 S 1/2 


766.49 


0.03173(7) 


0.4136(4) 





Theoretical results [35] from semiclassical model which does not resolve the doublet. 
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Table 5. Lorentzian half half-widths (in units of 10 21 MHz m 3 /atom) for alkali 
doublets pressure broadened by He perturbers. 



Element Transition T(K) Theory Experiment 



Li 




2 Pl/2 


- 2s 


2 S, to 
°l/2 


600 


0.3274 


0.302 a 


Li 




2 P 3 /2 


- 2s 


2 Si /o 
°l/2 


600 


0.3274 


0.302 a 


Li 


2p 


2 Pl/2,3/2 ~ 


2s 2 Si /2 


673 


0.3431 


0.33±0.17 b 


Na 


3p 


2 Pl/2 


- 3s 


2 s 1/2 


450 


0.2731 


0.260 ±0.003 c 
















0.242 ±0.016 d 


Na 


3p 


2 P 3 /2 


- 3s 


2 s 1/2 


450 


0.3096 


0.318 ±0.006 c 
















0.258 ±0.016 d 


Na 


3p 


2 Pl/2 


- 3s 


2 s 1/2 


480 


0.2807 


0.302 ±0.016 c 


Na 


3p 


2 P 3 /2 


- 3s 


2 s 1/2 


480 


0.3174 


0.349 ±0.032 e 


Na 


3p 


2 Pl/2,3/2 ~ 


3s 2 S 1/2 


2450 


0.567 


0.560 ±0.064s 
















0.726 ±0.064 h 


Na 


3d 


2 D 3/2 


-3p 


2 Pi/ 2 


470 


0.6290 


0.645 ±0.02 f 


K 


4p 


2 Pl/2 


- 4s 


2 s 1/2 


410 


0.2773 


0.247 a 


K 


4p 


2 P 3 /2 


- 4s 


2 Si/ 2 


410 


0.3808 


0.328 a 




e [35] 
f US 

s m 

h [42] 



Table 6. Lorentzian shifts (in units of 10 21 MHz m 3 /atom) for alkali doublets pressure 
broadened by He perturbers. 



Element 


Transition 






T(K) 


Theory 


Experiment 


Li 


2p 


2 Pi/ 2 - 


2s 


2 Sl/ 2 


630 


0.0023 


-0.045 a 


Li 


2p 


2 P 3 /2 - 


2s 


2 s 1/2 


630 


0.0011 


-0.018 a 


Na 


3p 


*l/2 


3s 


2 Sl/ 2 


450 


-0.0107 


0.048 ±0.014 b 


Na 


3p 


r 3/2 


3s 


2 Sl/ 2 


450 


-0.0139 


0.018 ±0.014 b 


Na 


3p 




3s 


2 Sl/2 


480 


-0.0105 


0.006 ± 0.009 c 


Na 


3p 


2 P , 
*3/2 ~ 


3s 


2 Si/ 2 


480 


-0.0134 


0.011 ±0.011 c 


Na 


3d 


2 D 3/2 - 


3p 


2 Pl/ 2 


470 


0.0352 


0.034 ± 0.003 d 


K 


4p 


2 Pl/2 " 


4s 


2 s 1/2 


410 


0.0988 


0.072 ±0.006 e 


K 


4p 


2 P 3 /2 - 


4s 


2 s 1/2 


410 


0.0250 


0.038 ±0.003 e 



a [13] 
b [37] 
c [H] 
d gD] 

e m 



